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ABSTRACT 



We investigate the effect of stellar activity and flares on short-term radial velocity measurements in the mid-M flare star CN Leo. 
Radial velocity variations are calculated from 181 UVES spectra obtained during three nights. We searched for spectral orders that 
contain very few atmospheric absorption lines and calibrated them against the telluric A-band from O2 in the Earth's atmosphere. 
One giant flare occurred during our observations, which has a very strong effect on radial velocity. The apparent radial velocity shift 
due to the flare is several hundred ms"' and clearly correlated with Ho- emission. Outside the flare, only spectral orders containing 
the most prominent emission lines of H, He, and Ca show a correlation to chromospheric activity together with a radial velocity jitter 
exceeding a few 10ms"'. We identify a number of spectral orders that are free of strong emission lines and show no flaring-related 
radial velocity jitter, although flares occurred as strong as 0.4 dex in normalized Ha luminosity. The mean radial velocity jitter due 
to moderate flaring is less than 10ms"'. Strong flares are easily recognized directly in the spectra and should be neglected for planet 
searches. 
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1. Introduction 

Radial velocity searches have led to a great many detections of 

■ extrasolar planets. Most of these planets are discovered around 
' stars of spectral types F, G, and K, although radial velocity am- 
plitudes are greater in stars of later spectral type and lower mass. 
The two reasons for the lack of planet detections around M 
dwarfs are a) the intrinsic faintness of M dwarfs render the de- 
tection of radial velocity variations very difficult, and b) many 
M dwarfs are very active, which is thought to be a problem for 
radial velocity detections, and it is sometimes thought to make 
M dwarfs less attractive in the quest for habitable planets. 

At the time of writing, 13 planets have been detected around 
i 9 M dwarf^. One of them, Gl 58 1 , is hosting two planets that po- 
' tentially fall into t he star's ha bitable zone (von Bloh et al. 2007; 
i ISelsis et alj|2007l) . Recently, iTarter et all 120071) has concluded 

■ that M dwarfs should be included in the quest for habitable 
' worlds and evidence of life, because they are so numerous and 

there are no obvious reasons why activity should destroy life on 
them. With the growing importance of M dwarfs, active or in- 
active, it is important to understand the systematic difficulties in 
detecting radial velocity variations. 

With the radial velocity technique, one detects planets by 
measuring radial velocity variations due to the orbital reflex mo- 
tion of the host star, which is typically on the order of a few 
m s~' depending on stellar and planetary mass. Radial velocities 
in general are determined by measuring the effective blue- or 
redshift of a spectrum relative to a reference spectrum. Such a 
shift in wavelength can either be imposed by a real velocity vari- 
ation or by a change in the appearance of the star's spectrum. 



* Based on observations collected at the European Southern 
Observatory, Paranal, Chile, 077.D-0011 
** Emmy Noether Fellow 
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The latter may stem from a cool spot on the surface of a star ef- 
fectively covering a part of the star's surface. Depending on the 
surface velocity on a rotating star, a spot can cau se virtual ra- 
dial ve locity variations well above the ms"' level. ^ Desort et al.l 
(I2OO7I) simulated the effect of starspots on the surface of F-, G- 
, and K-type stars rotating at moderate rotation velocities (up 
to V sin/ = 7kms"'). They find that starspots can lead to vir- 
tual radial velocity shifts of several ten ms"'; a spot covering 
a fraction of 1.07% of the surface of a G2 dwarf rotating at 
V sin / = 3 kms"' causes a radial velocity variation with an am- 
plitude of 60ms"'. 

The main effect of cool starspots is to remove (or dimin- 
ish) the flux coming from the part that is covered with a spot. 
Cool spots on the Sun are regions of strong magnetic fields that 
are several hundred to a thousand K cooler than their suiTound- 
ings. The higher atmospheric layers are heated by magnetic en- 
ergy so that in the vicinity of cool spots chromospheric lines 
like Ha become visible in emission. In very cool stars like M 
dwarfs, not much is known about the properties of st arspots 
and th en' e ffect on radial veloc ity measurements. Bo nfils et al.l 
(I2OO7I) and fPemorv et al.l (I2OO7I) report radial velocity variations 
due to starspots around the planet hosting M dwarfs Gl 674 
and Gl 436. The rotational periods are 35 d and 48 d, respec- 
tively, and semi-amplitudes ~ 5 ms"' were me asured. This 
is roug hly consistent with estimates of from Saar & Donahue! 
(II997I) . although their estimates are for a G2 dwarf with spot 
temperature zero, which certainly differs significantly from the 
situation in M dwarfs. 

In very cool stars, the contrast between chromospheric emis- 
sion lines and the quiet surface is much higher than in hotter stars 
so that observation of activity is easier in cooler stars. A forest of 
emission lines becomes visible during strong flares in M dwarfs 
(see e.g. Fuhrmeister et al. 2008). The influence of active regions 
on the spectra of M dwarfs, however, has not been studied in a 
systematic way. One part of the problem is that simulating chro- 
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Fig. 1. Normalized Ha luminosity during the three nights of observation. Time coverage of the three half-nights is shown in the plot. 
The three nights are shown individually in the insets, where the variability in Ha can be seen. A huge flare is clearly visible at the 
beginning of the first night. 



mospheric emission is much more difficult than simulating cool 
starspots. Nevertheless, it is expected that strong activity like 
flares have strong impact on the radial velocity measurements 
in M dwarfs. 

iKiirster et alj (2003) measured radial velocity variations in 
the inactive M dwarf Gl 699 (Barnard's star). They found an anti- 
correlation between radial velocity and the strength of the Ha 
line. In the spectrum of Gl 699, Ha is always seen in absorption. 
Even in the "flare" spectrum that Kiirster et al. show. Ha is not 
seen in emission but the absorption is weaker than in the other 
spectra (an eff'ect tha t may a ctually be due to lower temperature 
not higher, see, e.g.. lShort & Dovle 1998). Kiirster et al. found 
a significant anti-correlation between their Ha activity level and 
radial velocity. The small activity variations of Gl 699 led to ra- 
dial velocity variations on the order of 10ms However, the 
radial velocity measured in the "flare" spectrum of Kiirster et al. 
does not show as high a radial velocity shift as could be expected 
from the correlation seen at lower Ha level. 

The fact that cool spots on sun-like stars can cause large ra- 
dial velocity variations, and the result that the Ha line strength 
(anti-)correlates with radial velocities even in a very inactive M 
dwarf raise the question to what extent activity and flaring events 
influence the detectability of real radial velocity variations due 
to orbital motion. This is particularly important since in mid- 
and late-M dwarfs the fraction of active dwarfs rises dramati- 
cally, and the maj ority of mid- and late-M dwarfs show Ha in 
emission jWest et al..,2004) . 

2. Data 

In a multi-wavelength study of the very active nearby M6 dwarf 
CN Leo (Wolf 359, Gl 406), we obtained high resolution spec- 
tra with UVES at the Very Large Telescope, Paranal, Chile. The 
purpose of this data was to investigate the multi-wavelength be- 
haviour of flares. Here, we investigate this data for radial veloc- 
ity variations. The target always shows strong Ha emission, and 
during the three nights it presented many small flares and one 
very big one. 

We obtained 181 spectra during three half-nights in 2006 
( May 19/20, 21/22, and 23 /24). The data are described in detail 
in lFuhrmeister et akl (l2008l) where a spectral atlas of a spectrum 
taken during a huge flare is presented. The original purpose of 
the data was to investigate flares, not to search for radial veloc- 
ity varations. Hence, no special calibrations were done and the 
spectra were taken without the iodine absorption cell. We used 



a slit width of I" (R ~ 40000) and the red arm was centered at 
830 nm providing wavelength coverage from 640 nm to 1008 nm 
in 32 orders. 

Three half-nights were granted to the project. Because of 
bad weather conditions, no observations were taken during the 
first half of the second half-night. During the other two nights, 
CN Leo was observed more than four hours without interrup- 
tion. Exposure times varied between 100 s and 300 s depending 
on weather conditions. Several flares occurred during our obser- 
vations, among them a very strong one that boosted the observed 
Ha emission by more than an order of magnitude (it saturated 
the detector so we cannot measure the total Ha emission in that 
spectrum). The flare is ea sily visible in the w hole optical region, 
it is discussed in detail in lFuhrmeister et al.i (2008) . 

In Fig.[Tl normalized Ha activity is shown during the ob- 
servations. Measured Ha equivalent widths were converted to 
Lua/L\,oi using theoretical PHOENIX spectra as a reference (see 
iReiners & Basrill2007i) . Fig. [T] shows the temporal coverage of 
the data. In the insets, the three half-nights are shown individu- 
ally. The strong variability of CN Leo is obvious from this plot, 
and the huge flare during the beginning of the first night stands 
out as an extraordinary event. 

3. Radial velocities 

The most successful method detecting extra-solar planets so far 
is the radial velocity method. The hunt for radial velocity pertur- 
bations induced by planetary companions became possible when 
this technique became sophisticated enough to push the stability 
of radial velocity measurements down to several ten ms"', i.e. 
typically only 1/100 of a pixel and less. To reach such an accu- 
racy, wavelength calibration data is included in the target spec- 
trum itself, so that any instrumental effects are directly visible 
during the observation. One way to reach this is to obtain the 
spectrum of a calib ration lamp next to the target spectrum (the 
"ThAr method", see lPepe et al.ll2002h . A second method uses an 
absorption cell in front of the spectrograph. Light from the star 
passes the cell, which imprints a dense forest of absortion lines 
onto the spectrum (IButler et al.lll996l) . The crucial advantage of 
both methods is that every systematic effect in the wavelength 
solution of the spectrum is directly monitored during observa- 
tion (although the ThAr method requires much higher stability). 

Our spectra were not taken with special care regarding to 
wavelength stability. For calibration of UVES spectra, a ThAr 
spectrum is usually taken at the end of every night. These are the 
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Fig. 2. Relative radial velocities Avrad between the 181 spec- 
tra and an arbitrary reference frame (#24) versus local time at 
Paranal. Black crosses are data from the first night, blue trian- 
gles are from the second night, and green circles are spectra 
taken during night 3. Left: Radial velocities calculated in the tel- 
luric A-band; right: Av,ad calculated from order #69 that contains 
mainly the TiO IR band (see text). 



wavelength references that we use to calibrate our data. The typ- 
ical stability of UVES, however, is on the order of a pixel over 
the course of a night, i.e. deviations on the order of a kms ' 
can be expected. In order to reach an accuracy better than the 
absolute stability of the UVES spectrograph, we employ the ab- 
sorption lines from Earth's atmosphere. They provide a reliable 
reference frame on the order of the wind speed, w hich does not 
vary abruptly during a night dBalthasar et al.llT982h . None of the 
telluric lines were removed from the spectra prior to analysis. 

3.1. Relative radial velocity shifts 

To measure radial velocity shifts from our spectra, we first calcu- 
lated the cross correlation functions relative to a reference spec- 
trum. Spectrum number 24 obtained during the first night was 
used as reference spectrum; it exhibits relatively little Ha emis- 
sion. This was done for all spectral orders individually, i.e. we 
calculated 181x31 individual radial velocity diff'erences. In the 
following, these are labelled the relative radial velocity shifts 
(relative to an arbitrary zero point). They are not corrected for 
instrumental effects that may have occurred during the course of 
each night. 

Fig. |2] shows all 181 relative radial velocities Av,ad. Black 
crosses are data taken during the first night, blue triangles are 
taken during the second night, and green circles are observa- 
tions from the third night. In the left panel, Aviad from order #80 
(757-769 nm) is shown. This order covers the A-band due to 
atmospheric O2. Assuming no or very little variations in the tel- 
luric lines, Avrad from this order shows the absolute drift of the 
spectrograph during the nights. In the right panel of Fig.|2] we 
show radial velocity differences calculated from spectral order 
#69 that contains the IR TiO band and only very little telluric 
absorption. No barycentric correction was applied to the data 
shown in Fig.|2] 

Clearly, relative radial velocities of telluric lines and stellar 
lines follow the same gross pattern. During the first night (black 
crosses), the telluric and the stellar spectra show a strong red- 
ward shift within the first 15 minutes. More precisely, a blueshift 



occurs during the first three exposures followed by a steep red- 
ward jump. After that, a linear redward trend drifting approx- 
imately Ikms"' in four hours occurred. This trend, observed 
both in the atmospheric and the stellar lines, clearly is intro- 
duced by the UVES instrument itself and cannot be due to the 
motion of the star. The origin of the jump and the linear redward 
trend during the first night is unclear, but a thorough analysis of 
instrumental effects goes beyond the scope of this paper. During 
the second night (blue triangles) and the third night (green cir- 
cles), a shallow blueward trend is observed. At the beginning of 
the third night, a redward shift occurred that is similar to the one 
observed in the first night. The blueward drift during the sec- 
ond and third nights is on the order expected from changes due 
to variable atmospheric pressure inside the instrument (see next 
section). 

3.2. Differential radial velocity shifts 

Calibrating radial velocity shifts with telluric fines has several 
caveats that must be taken into account. The basic problem is 
that telluric lines are spread over almost the entire spectrum, but 
in particular in the rich spectra of M dwarfs they are not always 
easily discovered. The cross-correlation function of two spectra 
taken from the same star generally consists of two components: 
1. The correlation peak from stellar features, which follow the 
motion of the star and the drift due to barycentric motion. 2. 
The correlation peak of the telluric lines that are not subject to 
barycentric velocity drifts. If a spectral range contains both sets 
of features, the influence of the correlation peak due to telluric 
lines will affect the peak from stellar features. In the easiest case, 
telluric features dominate and the correlation peak no longer fol- 
lows the star's motion but rather the effects of telluric lines. In 
a more complicated situation, the varying distance between the 
stellar and the telluric peaks interfere so that the velocity mea- 
sured from the cross correlation function is a complex combina- 
tion of stellar and barycentric motion. 

One way to get rid of telluric features is to correct for them 
by using either a telluric standard spectrum or a synthetic set 
of telluric lines. Alternatively, one can mask out regions of tel- 
luric lines and use only regions that are virtually free of telluric 
contamination. As mentioned above, it is not always straightfor- 
ward to i dentify telluric feat ures in particular in the spectra of M 
dwarfs. Reiners et al. |(|2007) show a high resolution spectrum of 
a standard star in which telluric absorption at Paranal is visible in 
the wavelength region used for this work. This spectrum shows 
that telluric contamination is relatively weak outside the strong 
bands of O2 (A- and B-band at 760 nm and 690 nm, respectively) 
and the H2O features at 720 nm, 820 nm, and between 890 nm 
and 980 nm. This means that probably a large range of the spec- 
tral coverage of UVES can be used to measure the radial velocity 
shift of CN Leo. 

In this work, we chose neither to remove nor to mask tel- 
luric features. Instead, we computed the cross correlation func- 
tion for each order individually and search for orders that trace 
the barycentric velocities. We assume that the velocities calcu- 
lated from these orders follow the apparent radial velocity of the 
target. 

With the relative radial velocity shift with respect to one ref- 
erence frame, we can now calculate the differential radial ve- 
locity shift Av'rad.stai- " Av'rad.teiiuric - Avs - Avt, i.e. the rcsidual 
between the relative telluric shift from order #80 (left panel of 
Fig.|2]i and the relative stellar shift (one order shown in the right 
panel of Fig.|2]i. The differential radial velocities for all orders 
are shown in Fig.|3] We group the orders into "stellar" and "tel- 
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Fig. 3. Differential radial velocities without barycentric corrections for all spectral orders individually. Black crosses: first night; 
blue triangles: second night; green circles: third night. Top three panels show orders that resemble the pattern of barycentric motion. 
These orders are compiled in Table[T] The barycentric motion is shown in the upper right plot. Lower two panels show spectral 
orders not following barycentric motion. These orders are affected by telluric lines or other contaminants, they are compiled in 
Tabled 



luric" orders based on whether or not they trace the barycentric 
motion. The three top panels show "stellar" orders following the 
apparent velocity introduced by the motion of Earth, which it- 
self is shown in the upper right panel of Fig. [3] These orders are 
compiled with their wavelength coverage in Table[T] The am- 
plitude of the drift introduced by the motion of Earth is on the 
order of 500ms ' per 4h and clearly dominates the radial ve- 
locity signal measured during each night in the 20 orders shown 
here. The radial velocity differences between the three individual 
nights (sets of black crosses, blue triangles, and green circles) do 
not strictly resemble the pattern expected from barycentric mo- 
tion. This is not suprising because the absolute stability of the 



spectrograph is not better than a few hundred ms ' from night 
to night. The differences between the first and third night are in 
general closer to the expectations than the differences between 
the second and the other two nights. The differences vary be- 
tween different orders. We calculate the median of the differen- 
tial radial velocity shifts from the first and second nights and 
plot their differences in Fig.|4] as a function of spectral order 
The offset between spectra taken during the first and the second 
night varies between ~ -500 and 300 m s ' and shows relatively 
smooth variation with wavelength, which indicates that the off- 
set is introduced by inaccuracies of the dispersion relation. Thus, 
UVES data calibrated with only one ThAr frame per night can 
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Table 1. Spectral orders of UVES that contain mostly stellar ab- 
sorption features. Orders with relevant emission lines or telluric 
absorption lines are excluded. 
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Fig. 4. Offset between the median differential radial velocities 
during night #2 and night #1 as a function of spectral order 
(difference between blue triangles and black crosses in Fig. [3). 
Filled circles are from orders tracing the star's motion, open 
stars show telluric orders (top and bottom panel in Fig.|3), re- 
spectively. 



Table 2. Spectral orders of UVES not used for radial velocity 
calculations. Telluric absorption features or broad stellar features 
are indicated. 



Order # 


Wavelength coverage [nm] 


feature 


62 


976-991 


H2O 


63 


960 - 976 


H2O 


64 


945 - 960 


H2O 


65 


931 -945 


H,0 


66 


917-931 


H,0 


67 


903-917 


H2O 


68 


890 - 903 


H2O 


70 


865 - 877 


Ca (stellar) 


75 


808-818 


HoO 


80 


758 - 767 


02 


86 


705-713 


He, TiO (stellar) 


89 


681 -689 


O2 


93 


652-659 


Ha (stellar) 



clearly not be used to search for radial velocity variations on the 
order of 100m s ' over several hours. On the other hand, short- 
time radial velocity variations can be traced as will be shown in 
the following. 

The two bottom panels of Fig. [3] show differential radial ve- 
locities in the spectral orders that do not follow the drift expected 
from barycentric motion. These are the orders that are either af- 
fected by large scatter caused by broad features not suitable for 
precision velocity measurements (orders #86 and #93), or that 
are affected by telluric contamination. Order #70 was taken out 
because it follows the weaker flare at the beginning of night 3 
(before 20:00), which is probably caused by the Ca emission 
line. The spectral orders shown in the bottom two panels are not 
used for radial velocity calculations, they are given in Table|2] 
together with the source of contamination. 

A comparison between the spectral orders of Table[T] and 
known telluric absorption features shows that some spectral re- 
gions, e.g., at 690 nm and at 720 nm contain telluric absorption. 



0.3 r 




-0.1 r 



20;00 22:00 00:00 

Local Time 

Fig. 5. Differential radial velocities Avs - Avt for the three nights 
versus local time. Symbols are the same as in Fig.|2l A second 
order polynomial fit for the three nights is overplotted in red. 
The giant flare at the beginning of the first night was excluded 
from the fit. 



These orders are likely to be affected by telluric absorption, too, 
which can be a problem for radial velocity work. At least they 
should be used only with great care. In our case. Fig. |3] shows 
that the stellar signal is dominating the radial velocity measure- 
ment (e.g., in orders 84 and 88), and we chose to include them 
in the set of "stellar" orders. Furthermore, a comp arison with the 
telluric standard spectrum in lReiners et al.l ( l2007h . taken with the 
same instrumental setup, shows that the telluric features identi- 
fied in Table|2] represent the strong features seen from Paranal. 
No other significant absorption features (on the order of 20 % 
depth) appears in this atlas. The relevant parts of the telluric 
bands at 690 nm and 720 nm, for example, have less than 15 % 
depth, which is much weaker than the depth of stellar features. 
We note that classifying more orders that are suspect to telluric 
contamination as "telluric" has no significant influence on the 
results presented in the following. 
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Fig. 6. Differential velocities plotted versus atmospheric pres- 
sure during the three nights of observation. Symbols are the same 
as in Fig|2] 



Our next step is to apply the barycentric correction and to 
calculate the mean differential radial velocity shift from the 
"stellar" orders (top three rows of Fig.[3]and Table[Tll- The mean 
of the differential velocity shifts from the 20 orders given in 
Table[T]is shown in Fig.|5] 

The differential radial velocity shifts from each of the three 
nights all lie on smooth curves. Differential velocity shifts from 
the second and the third night show very similar absolute val- 
ues. Those from the first night, however, fall about 200ms ' 
below the results of the other two nights. This discrepancy is 
probably connected to the strong drift in the relative radial ve- 
locities during the first night. We checked temperatures during 
the three nights and found that relatively large temperature vari- 
ations occurred during the first night. We can speculate that this 
may be connected to the systematic offset between the disper- 
sion relations during that night and the other two nights (Fig.|5]). 
However, the temperature variations cannot explain the large 
drift in relative velocities seen in Fig. |2] (the temperature effect 
is expected to be much smaller, see Kaufer et al. 2007). 

The smooth variation of differential radial velocity shifts ob- 
served during each of the three nights are all following a very 
similar pattern. It is not the main objective of this paper to inves- 
tigate the nature of this smooth pattern, because we are mainly 
interested in short term variations connected to stellar flares that 
we can trace in Ha. Nevertheless, we can speculate about the 
reason of the smooth changes seen in each of the three nights. 
The striking similarity of the nightly variations suggests a rea- 
son connected to the nightly course of observations. One possi- 
bility is that some interference between telluric lines and stellar 
lines is left in the orders that we selected above. As explained, 
this could introduce changes of radial velocity drift with chang- 
ing barycentric motion. However, the barycentric motion is al- 
most monotonous during each night (upper right box in Fig.|3ll. 
The nightly velocity variations consistently show a maximum 
around 23:00 h each night, which makes a connection to telluric 
lines improbable (although we cannot completely exclude this 
option). 

A second possibility to explain the intra-night pattern are 
variations of atmospheric pressure inside the UVES instrument. 
Pressure variations can alter the dispersion relation and hence 
affect the differences between wavelength shifts observed at dif- 
ferent spectral orders. We checked the behaviour of atmospheric 
pressure during the three nights; they indeed follow a pattern 
very similar to the curves in Fig.|5] In Fig.|6] the differential ve- 
locity shift is plotted as a function of pressure. For each night, 
the curves follow a relative well defined line with a similar de- 
pendence on pressure while the absolute pressure values are very 



different. The amplitude of pressure varia tions is consistent w ith 
the expectations from the UVES manual (iKaufer et al.ll2007l) . 



3.3. Could the intra-night drift be reai? 

A viable explanation for the smooth trend in differential radial 
velocities in Fig.|5]is a change of dispersion coupled to atmo- 
spheric pressure inside the instrument. However, this explana- 
tion cannot be tested, and the trend may as well be real. If 
we assume that the offset between the differential radial veloc- 
ities, measured during the first and the other two nights, is in- 
strumental but the drift comes from real radial velocity varia- 
tions, the "detected" period would be about 2 days or an inte- 
ger fraction of that. A potential explanation for such a period in 
the flare star CN Leo are corotating spots producing radial ve- 
locity variations with a semi-amplitude of Ks ~ 100 ms"'. As 
in higher mass stars, spots may produce ra dial velocity varia- 
tions on the order of several 10 to 100 m s~' jBonfils et 31^20071 : 
iDemorv et al.ll2007l;lDesort et alj|2007h . The rotational period of 
CN Leo is not kno wn, but it s surface rotat ion ve locity is about 
V sin/ = 3kms"' dReiners & Bas ri 2007). lEacvl ( [197 7) gives a 
value of logR/RQ = -0.84 {R = Q.URq) for flie radius of 
CN Leo. Thus, the (projected) rotational period of CN Leo is 
P/ sin/ «i 2d! This means that rotational modulation as a rea- 
son for the radial velocity changes seen during each of the three 
nights is not inconsistent with the expected rotational period of 
CN Leo. However, we have tested the periodicity of Ha and 
other emission lines and found no sign of a rotational period of 
2 d. This may imply that the rotational variation is buried under 
the strong and irregular flaring activity so that CN Leo's rota- 
tion cannot be detected in emission lines, but it may also that the 
rotational period of CN Leo is very different from 2 d. 

A third way to generate radial velocity variations as observed 
in our data is an object in orbit around CN Leo. The semi- 
amplitude caused by the hypothetical planet would be around 
1 00 m s ' , which could be explained by a planet with a projected 
mass of Msin/ ~ 0.2 Mj at a distance of 0.015 AU. For now, 
this hypothesis cannot be falsified with the available data, but 
seems to be extremely unlikely given the close match to pressure 
changes in the instrument and the coincidence with the estimated 
rotational period of CN Leo. 



4. Results 

The main focus of this work are the short-term variations of the 
apparent radial velocity introduced by stellar flares. These flares 
have timescales of the order of about half an hour 

We removed the intra-night trend from the differential radial 
velocity shifts in Fig.|5]by applying a second order polynomial 
fit (exluding the large flare during the first night). The resulting 
final radial velocities are discussed in the following. 



4. 1 . Radiai veiocity jitter 

Fig.|2]displays the rms scatter for all UVES orders that are cov- 
ered by our setup. Again, the six spectra taken during the huge 
flare were neglected. Only nights 1 and 3 are shown because 
much less data is available for the second night. Spectral or- 
ders affected by telluric lines (Table|2]i are marked with a star 
Symbols as in the figures above, i.e. black crosses for night 1 
and green circles for night 3, are used for the other orders cover- 
ing mainly stellar features (Table[T]i. These orders are connected 
by black and green lines for the first and third night, respectively. 
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Fig. 8. Radial velocity jitter for the three nights (from top to bottom, symbols as above) together with normalized Ha luminosity 
(grey circles). The polynomials shown in Fig.|5]are subtracted from the differential radial velocities. 
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Fig. 7. Scatter of differential radial velocities around the trend 
fitted in Fig.|5] Only nights 1 (black crosses) and 3 (green cir- 
cles) are shown. Stars indicate orders that predominantly contain 
telluric features, arrows mark the orders used to compute stellar 
radial velocities (these are the orders given in Table[TJ. 



Interestingly, the rms scatter is decreasing among the lowest or- 
ders (longest wavelengths) from values around 30m/s in order 



#61 down to values around 10 m/s in order #69. The reddest or- 
der is virtually free of telluric features but contains many strong 
stellar lines due to molecular FeH. The next 7 orders predomi- 
nantly contain absorption due to atmospheric H2O. The trend of 
higher scatter in the reddest orders is probably caused by the 
quality of the wavelength solution, which becomes worse to- 
wards the near infrared because the ThAr calibration lamp does 
not provide as many lines as in the bluer spectral orders. 

Most of the orders have rms values between 10 and 20 m s"' . 
Several orders show relatively high rms values. This "scatter" 
is correlated with Ho' emission (see next section) and is caused 
by chromospheric emission lines contained in these orders. The 
most prominent lines that affect radial velocity measurements 
are the lines of the Ca infrared triplet at 849.8 nm, 854.2 nm, and 
866.2 nm (orders #72, #71, and #70), two He I Hnes at 706.5 nm 
(order #86) and 667.8 nm (orders #91 and #92), and Ha (order 
#93). Order #86 also contains the strong TiO y-band. This very 
strong feature is a dense forest of blended lines generating wide 
correlation functions, which are part of the reason for the rela- 
tively large scatter in radial velocity. 

In Fig. [8] the mean differential radial velocity shifts calcu- 
lated from all orders in Table[T] are shown after substracting 
the smooth trend (Fig.|5]l. The first night is plotted in the top 
panel, the second in the middle, and the third night in the bottom 
panel. In the same figure, the values of normalized Ha emission 
logLiia/Lboi are plotted as grey circles. 
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Fig. 9. Correlation coefficients between the radial velocity shifts 
and normalized Ha luminosities. Black crosses (+) and green 
circles show correlation coefficients for each order containing 
mainly stellar features during the first and third nights, respec- 
tively (Table[T]i. Stars show orders of Table|2l As red crosses 
(x), the correlation coefficients calculated for the first night are 
shown including the spectra taken during the huge flare. 



4.1 .1 . The huge flare 

In Fig. [8] the strong flare that occurred at the beginning of the 
first night is clearly visible as a spike in radia l velocity. Spectra 
of this huge event are discussed in detail in iFuhrmeister et alj 
(12008V such a flare is obvious in the spectra of M dwarfs. In 
the spectrum at flare peak, the effect on radial velocity is as 
large as several hundred m s"' and it shows the opposite sign 
(Avs - Avt - -660ms"', not shown in Fig.©. The exact value 
of the radial velocity shift during flare peak has nothing to do 
with the star's radial velocity (it may reflect the geometry and 
dynamics of the flare). We have not included this point in Fig. [8] 
but we note that its radial velocity is an order of magnitude larger 
than in all other spectra. After the flare, the radial velocity shift 
diminishes together with Ha luminosity. 

In a search for radial velocity variations, the spectra affected 
by a flare as huge as this one are easily distinguished from 
"quiet" spectra. In a planet search program, they should not be 
used to detect radial velocity variations. 

4.1 .2. Jitter outside the huge flare 

The radial velocities outside the huge flare show no obvious 
trend or period. The rms scatter for the three nights is 8.2 m s ' 
for the first night, 8.4ms ' for the second night, and 9.7ms ' 
for the third night. The overall jitter for all three nights taken 
together is 9.0ms The meaning of this result is twofold: a) 
For the purpose of this investigation, the accuracy to which dif- 
ferential radial velocity variations are determined is on the order 
of 10 ms"'; b) outside the giant flare, the still relatively strong 
flaring activity 0.4 dex in log Letr/Lboi) on CN Leo has no 
effect on radial velocity measurements in excess of 10ms "' as 
long as spectral orders are used that avoid strong chromospheric 
emission lines. 



4.2. Correlation with activity 

In the calculation of the rms jitter above, we excluded the spectra 
taken during the huge flare. This flare event does show influence 
on the radial velocity curve, but flares of this strength pose no 
danger to radial velocity surveys because they are obvious in 
the spectra; they must not be used for planet searches. However, 
the data on Ha emission shows that a number of smaller flaring 
events occurred during the three nights. Are the smaller flares, 
which are not easily detectable in the spectra and will contribute 
to radial velocity measurements, affecting the radial velocity 
measurements? From a first glance, the curves of radial velocity 
and Ha luminosity in Fig.[8]show no striking similarity. To check 
the correlation between radial velocity and Ha activity, we com- 
puted the linear Pearson correlation coefficients (IDL procedure 
correlate) between normalized Ha luminosity and differential 
radial velocity. 

The results for each spectral order are shown in Fig.|9] 
Symbols are the same as in Fig.|2] Data for the first night are 
computed without the six spectra taken during the strong flare. In 
most of the spectral orders, the absolute value of the correlation 
coefficient is small (< 0.3) indicating no or very weak correlation 
between Ha emission and radial velocity. Thus, flaring shows 
no obvious effect on the radial velocities as long as Av is mea- 
sured in spectral orders free of chromospheric emission lines. 
On the other hand, spectral orders containing prominent emis- 
sion lines, namely orders #70, #72, and #93, exhibit rather high 
(absolute) correlation coefficients. Here, activity strongly affects 
radial velocity measurements through chromospheric emission 
in the Ca, He, and Ha lines, although small flares cannot imme- 
diately be detected in the spectra themselves. These orders must 
to be omitted when searching for radial velocity variations in 
active M dwarfs. 

In Fig.|9] we also show the correlation coefficients calcu- 
lated from the entire first night, i.e., including the huge flare 
(plotted in red). As mentioned above, this flare substantially af- 
fected the whole spectrum including orders far away from Ha 
(see Fuhrmeister et al.. 2008) . The red crosses in Fig. |9] demon- 
strate that almost all spectral orders - including the ones con- 
taining mostly telluric lines - now show strong correlation with 
Ha emission. Thus, the radial velocity shift measured during this 
flare clearly is a systematic effect and not due to motions of the 
flare plasma. Spectra taken during flare events similar to this one 
therefore carry little information about the real radial velocity 
of the star (on the order of a few hundred ms"'). Interestingly, 
the two orders containing the Ha line show an anti-correlation 
with Ha emission. An anti-correlation between radial velocity 
and Ha luminosity was also measured at h igher accuracy in a 
much less active star bv lKiirster et al.l (l2003l) . 

5. Summary 

We have used UVES observations of the very active M dwarf 
CN Leo during three nights investigating the effect of flares 
on short-term radial velocity determinations. Radial velocities 
were measured relative to telluric absorption bands because no 
bracketing ThAr exposures or calibrations through an iodine cell 
were available. With this method, we could remove strong abso- 
lute shifts that occurred during the nights. Remaining intra-night 
drifts may be explained by pressure variations changing the dis- 
persion of the spectrograph or by corotating spots on the surface 
of CN Leo. These smooth trends were neglected in our analysis. 

During two nights, data were taken continuously for more 
than four hours. This allowed to investigate short-term variations 
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of radial velocity. Many flaring events were observed during that 
time. The timescale of flares is less than an hour and several 
flares are covered by our data. 

At the beginning of the first night, a huge flare occurred 
that dramatically changed the appearance of the whole spectrum. 
This flare had strong effect on the measured radial velocities in 
all spectral orders. Even spectral orders that predominantly con- 
tain telluric lines show radial velocity shifts on the order of a few 
hundred ms ^ This dramatic kind of flaring events, however, is 
not a general problem for radial velocity surveys, because such 
an event is easily identified in the spectrum. In this case, the 
flare limits the radial velocity stability to ~ 500 m s"' in the peak 
spectrum, and to ~ 50ms ' directly after the peak. Such events 
must be removed from the analysis. 

Outside the huge flare, spectral orders containing prominent 
emission lines show large scatter. Radial velocities measured 
in the order containing the Ho- line are anti-correlated to Ho- 
emission strength. This was also repor ted in the very inactive 
M dwarf Gl 699 bv lKursteretaD(l2003h . 

In the orders free of prominent emission lines, the radial ve- 
locity jitter is below 10ms Outside the huge flare, CN Leo 
still showed flaring activity that enhanced the normalized Ho- 
flux by ^ 0.4 dex. Despite this rather strong activity, no correla- 
tion was found between activity and radial velocity. 

Even in one of the most active known M dwarfs, flares pro- 
duced during average stellar activity have no effect on high pre- 
cision radial velocities on the level of 10ms ' if one chooses 
spectral orders free from emission lines. This means that the in- 
fluence of flares - the most obvious effect of stellar activity on 
the spectra of M dwarfs - is negligible in planet searches at least 
down to several m s"' . Even in the presence of flaring, it may be 
possible to detect an earth-Uke planet inside the habitable zone 
of a late M dwarf. 
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